Electron quantum dynamics in closed and open potentials at high magnetic fields: 
Quantization and lifetime effects unified by semicoherent states 
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We have developed a Green's function formalism based on the use of an overcomplete semicoherent 
basis of vortex states specially devoted to the study of the Hamiltonian quantum dynamics of 
electrons at high magnetic fields and in an arbitrary potential landscape smooth on the scale of 
the magnetic length. This formalism is used here to derive the exact Green's function for an 
arbitrary quadratic potential in the special limit where Landau level mixing becomes negligible. This 
solution remarkably embraces under a unified form the cases of confining and unconfining quadratic 
potentials. This property results from the fact that the overcomplete vortex representation provides 
a more general type of spectral decomposition of the Hamiltonian operator than usually considered. 
Whereas confining potentials are naturally characterized by quantization eff'ects, lifetime effects 
emerge instead in the case of saddle-point potentials. Our derivation proves that the appearance of 
lifetimes has for origin the instability of the dynamics due to quantum tunneling at saddle points 
of the potential landscape. In fact, the overcompleteness of the vortex representation reveals an 
intrinsic microscopic irreversibility of the states synonymous with a spontaneous breaking of the 
time symmetry exhibited by the Hamiltonian dynamics. 
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I. INTRODUCTION 

The integer quantum Hall effect, with its remark- 
able transport properties^ offers perhaps the simplest 
route to understand the complex dynamics of electrons 
taking place in strongly inhomogeneous nanostructures. 
Indeed, the presence of a strong perpendicular mag- 
netic field in two dimensions brings the classical mo- 
tion close to integrability, with slow drift trajectories 
superimposed to faster cyclotron orbits, leaving hope 
that the quantum dynamics displays similar and simple 
structures i^i^ii^i^i^i^ii^ii^ Further simplification is brought 
by the fact that electron-electron interactions can be 
taken into account in the integer quantum Hall regime 
at the single-particle level,— so that the calculation of 
equilibrium properties, such as the local electron density 
and the distribution of permanent currents throughout 
the sample, can be carried out from a one-particle ran- 
dom Schrodinger equation. Nevertheless, the precise mi- 
croscopic resolution of this equilibrium one-body problem 
still lacks a concrete analytical formalism and progress 
in this direction would be useful toward a microscopic 
description that underlies more complex noncquilibrium 
phenomena. 

The main difficulty in the resolution of the disordered 
Schrodinger equation resides in the complexity of the po- 
tential landscape that competes with the kinetic energy 
in a nontrivial manner. In fact, it is worth mentioning 
that the standard procedurei^ to deal with a random po- 



tential, which consists in averaging over impurity config- 
urations, is already physically questionable at high mag- 
netic fields, at least at the microscopic level. Indeed, this 
theoretical route is usually justified by the physical as- 
sumption of randomness after successive collision events. 
However, instead of the chaotic exploration of the disor- 
dered landscape, the electronic classical motion becomes 
relatively regular for a smooth disorder potential at high 
magnetic fields. At the technical level, this difficulty was 
pointed out in the standard quantum-mechanical dia- 
grammatic perturbative method, which leads to unsolved 
complications for a smoothly disordered potential: more 
and more classes of diagrams must be incorporated in 
the calculation as the strength of the magnetic field is 
increasedJ^ This somehow indicates that the perturba- 
tive technique is unadapted for the high magnetic-field 
regime when the bending of the electronic trajectories 
becomes too important. 

We note already that if one does not resort to any aver- 
aging over disorder, one faces not only a technical prob- 
lem but also a more fundamental physical one, namely, 
the question of the microscopic origin of irreversibility 
and intrinsic dissipation from Hamiltonian quantum dy- 
namics, an essential aspect for the calculation of trans- 
port properties. Indeed, the standard impurity averaging 
procedure, valid in the limit of low magnetic fields, in- 
troduces an effective description of the formalism. This 
allows to get, in addition to the energy spectrum, life- 
time effects, which are inaccessible in a purely Hamil- 
tonian formalism limited to the Hilbcrt space of square 
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integrable functions. In absence of averaging, it is thus 
necessary to clarify the possible relation between the 
"complexity" of the potential landscape and the issue on 
the microscopic origin of irreversible processes. In other 
words, we are confronted with the controversial ques- 
tion whether irreversibility results from supplementary 
approximations to the fundamental quantum-mechanical 
laws (which arc strictly time reversible) or is subtly hid- 
den in the usual formulation of quantum dynamics. 

A well-known starting point to capture the regime of 
quasi-regular dynamics in high perpendicular magnetic 
fields is to implement in full exact quantum-mechanical 
terms the fast cyclotron motion of circling electrons re- 
sulting from the Lorcntz force, which gives rise to the 
quantization of the kinetic energy into discrete Landau 
levels. The second relevant degree of freedom then cor- 
responds to the slow guiding center of motion, whose 
dynamics is dictated by the smooth disordered potential 
landscape. The essential role of the potential landscape, 
to be captured precisely in the microscopic quantum- 
mechanical calculations, is to lift the huge degeneracy of 
the Landau levels. It is already important to note that 
the magnetic field B enters into the quantum-mechanical 
problem only via two different quantities: the cyclotron 
pulsation ujc = |e|i?/TO*c and the magnetic length Ib = 
y''hc/\e\B. While uJc is a material-dependent parame- 
ter via the effective mass m*, Ig can be regarded as a 
more fundamental quantity since it involves only phys- 
ical constants such as Planck's constant H, the speed 
of light c, and the absolute value of electric charge |e|. 
The cyclotron pulsation Uc determines a characteristic 
frequency for the circular motion which actually already 
enters into the problem at the classical level. In contrast, 
the magnetic length Ib is purely a quantum-mechanical 
quantity, especially characterizing the spatial extent of 
the wave functions. In the popular operatorial language 
of quantum mechanics, these two relevant degrees of free- 
dom are introduced by decomposing the electronic coor- 
dinate f = R -I- r) into a relative position t) = v x z/wc 
linked to quasicircular cyclotron orbits (v is the velocity 
operator), and a guiding center position R = {X,Y). 
These quantum variables obey the commutation rules 
['Sajj'yy] = —ihhJc/m* and = iPg. In close analogy 

with the canonical quantization rule between the position 
and the momentum, it is seen that, for the slow drift mo- 
tion of the guiding center, plays the role of an effective 
magnetic- field-dependent Planck's constant. 

The condition of a high magnetic field can thus be im- 
posed cither by expressing that Ib is the smallest length 
scale in the problem, i.e., by taking Ib — > 0, or by consid- 
ering that UJc is the biggest frequency scale in the prob- 
lem, i.e., by taking uj^ — > oo. These two limits can in 
principle be taken separately or simultaneously and ac- 
tually yield different physical situations which have been 
discussed in the literature.i^'^i'^'iiii^'i^iii'i^'i^ For in- 
stance, a popular approach corresponding to the first 
limit Ib — * is to treat classically the slow guiding center 
motion while the fast cyclotron motion is kept quantum 



mcchanicalj^i^i^iSiii^ This case leads to great simplifica- 
tions in the theoretical treatment, since the guiding cen- 
ter coordinates then commute and can be described en- 
tirely in classical terms. In this limit, the guiding center 
motion is restricted to equipotcntial lines and the en- 
ergy spectrum is characterized by a continuous potential 
energy on top of discrete Landau levels. Another stan- 
dard approximation corresponding to the second limit 
ciJc ^ oo is to neglect Landau level mixing and to study 
the Hamiltonian quantum dynamics at finite Ib projected 
onto a single Landau leveL '^i-'^^i-'^^i-'^'^i-'^^'-'^^ This most tricky 
regime implies to work in a fully quantum-mechanical 
formalism taking into account rigorously the noncommu- 
tativity of the guiding center coordinates. An interest- 
ing aspect that was not fully cleared with both types 
of approaches lies in how to capture the transition from 
quantum to classical, with the classical features emerging 
possibly from microscopic decoherence processes. 

In order to study the interplay of Landau levels quanti- 
zation and a smooth disordered potential in a controlled 
fashion, we have devclopcd^Si^ in recent years a spe- 
cially devoted Green's function formalism based on the 
use of a semicoherent overcomplete set of states |to,R) 
labeled by a continuous quantum number R, related to 
the classical guiding center coordinates, and an integer 
TO, associated to the discrete Landau levels. This family 
of states were named vortex states^S due to the vortex- 
like phase singularity of the associated wave functions 
(r|TO, R) at the electronic position r = R. Because the 
vortex states encode no preferred symmetry, they allow 
a great adaptability to the local spatial variations of the 
random potential. More precisely, our approach consists 
in mapping the quantum equation of motion obeyed by 
the Green's function (the so-called Dyson equation) to 
this vortex representation, which then rigorously extends 
to quantum mechanics the classical guiding center pic- 
ture. The main essential difference with the guiding cen- 
ter treatment is that our method keeps the full quantum- 
mechanical noncommutativity of the guiding center co- 
ordinates through the overcompleteness property of the 
basis of states. We have shown^l that, within the vortex 
representation, Dyson equation can be easily and system- 
atically diagonalizcd order by order in powers of the mag- 
netic length Ib- Quantum observables are then obtained 
by returning to electronic representation from the vortex 
Green's functions, so that the semiclassical limit ^ 
as well as its systematic corrections is naturally obtained 
with our approach^i Moreover, the vortex representa- 
tion allows one to classify and include in a systematic 
and straightforward way the Landau level mixing pro- 
cesses in the calculations. 

The results to be developed here aim at extending the 
work initiated in Refs. I20ll2ll . where a systematic and 
closed form expression for the solution of Dyson equation 
in a smooth arbitrary potential was already obtained, 
under the form of a series expansion classified order by 
order in powers oi Ib- The further and important step 
made in the present paper is to include to all orders the 
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contributions from first and second spatial derivatives of 
the potential, in loose analogy to the resummation of 
leading classes of Feynman diagrams in standard pertur- 
bation theory. For simplification, we will consider the 
mathematical limit where Landau level mixing can be 
considered negligible and our solution will encompass all 
cases of quadratic potentials in that limit. A further mo- 
tivation for a resummation of the gradient expansion is 
that the series obtained in powers of Ib may not con- 
verge in general, since the semiclassical guiding center 
limit — > is expected to be singular, similar to the 
case of the more standard, fully semiclassical limit ft ^ 0. 
That the small Ib expansion is indeed singular will be il- 
lustrated by the asymptotic character of guiding center 
semiclassical results: physical aspects related to an ex- 
act quantum treatment, such as quantization of energy 
levels or lifetime effects, can not be approximated by a 
finite expansion. We shall confirm this feature by com- 
paring our exact quantum solution to various approxi- 
mation schemes related to several improvements of the 
semiclassical guiding center method. 

Before diving into the heart of the technique, we want 
to mention that, although our method obviously shares 
on certain aspects some similarities with theories al- 
ready existing in the literature, important differences 
with these prior works can be emphasized. First, we 
would like to stress that our methodology is based on 
the exclusive use of Green's functions, not wave func- 
tions, in contrast to the theory pioneered by Girvin and 
Jacbi (see also Refs. [l6lll8lfl9[ ) where a one-dimensional 
(ID) Schrodinger's equation for the electron dynamics 
projected onto a single Landau level (valid in the limit 
Wc ^ oo) was derived for finite Ib- This point could 
appear superfluous at the first glance, since it is possi- 
ble to get Green's functions from the knowledge of wave 
functions. However, the use of an overcomplete represen- 
tation with nonorthogonal states to solve the dynamical 
equations of motion necessarily forces us to give up the 
wave functions picture and work in a Green's function 
formalism of partially coherent states. Furthermore, it is 
worth noting that the Hilbcrt space of square intcgrable 
wave functions is usually well suited for closed intcgrable 
systems, but turns out to be totally inadequate in sit- 
uations presenting scattering processes in open systems 
(case of a saddle-point potential for instance) , where one 
must appeal to another formalism, the scattering states 
picture. We shall show that the use of an overcomplete 
representation of coherent states allows one to get and 
treat quantization and lifetime effects on an equal foot- 
ing in the resolution of Dyson equation. Moreover, the 
appearance of lifetimes in the energy spectrum coincides 
with the impossibility to describe the solution of the 
Dyson equation in terms of a countable set of states, 
thus proving the relevance of an overcomplete^^ repre- 
sentation in such a situation. Second, we note that sev- 
eral authors2ii^ii^iii actually attempted to build theories 
based on the use of vortex states within the path-integral 
formalism, which, however, seemed to suffer from techni- 



cal difficulties that were not elucidated i^ '^^i"'^^ In contrast, 
our theory is not tainted with the specific mathemati- 
cal ambiguities which can often be encountered with the 
path-integral technique. 

The derivation of an exact solution for the Green's 
function at large cyclotron energy yet finite magnetic 
length, embracing all possible cases of quadratic poten- 
tials, constitutes the main mathematical result of this 
paper. Besides capturing exactly the tunneling processes 
in the case of a saddle-point potential, it has the virtue 
of pointing out clearly the physical microscopic mech- 
anism responsible for the appearance of lifetimes in the 
spectral decomposition of the Hamiltonian. We therefore 
hope that it will also help to clarify the debate about the 
physical roots of time irreversibility and the ubiquitous 
emergence of a classical character from quantum mechan- 
ics. An important point we will also demonstrate is that 
the derived solution provides a controlled approximation 
at finite temperature for all equilibrium local observables 
in the case of an arbitrary potential that is smooth on 
the scale of the magnetic length. This result, based on 
the fact that the local Green's function at high magnetic 
fields displays a hierarchy of energy scales controlled by 
successive spatial derivatives of the potential, can be used 
to write down an expression for the local density of states 
which may be useful in the context of recent scanning 
tunneling spectroscopy measurements 3^ A short report 
of parts of this work has been published in Ref. [2^ . 

The paper is organized as follows. In Sec. |ll]we present 
the vortex Green's function formalism and derive the gen- 
eral form of Dyson equation in the vortex representation. 
In Sec. mil Dyson equation is exactly solved for the two 
particular cases of an arbitrary ID potential and an ar- 
bitrary two-dimensional (2D) quadratic potential. The 
obtained solutions are then exploited in Sec. |lV]to derive 
a general expression for the local density of states. Fi- 
nally, we discuss in Sec. |V]thc importance of considering 
an overcomplete representation in the present problem 
and its physical implications for the issue of time irre- 
versibility. A small conclusion closes the paper. Some 
extra technical details bringing complements of informa- 
tion for the calculations are given in several appendixes. 

II. DYSON EQUATION IN VORTEX SPACE 

A. Hamiltonian and projection onto the vortex 
representation 

We consider the single-particle Hamiltonian for an elec- 
tron of charge e = — |e| confined to a two-dimensional 
(x, y) plane in the presence of both a perpendicular mag- 
netic field B and an arbitrary potential energy ^(r), 

H=-^ (-thVr - -A(r)) ' + V{r), (1) 
2m* V c / 

with the vector potential A defined by V x A = B = Bz, 
and m* the electron effective mass (here r = (x, y) is the 
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position of the electron in the plane) . 

For V = 0, the energy spectrum is quantized into Lan- 
dau levels Em = [m + l/2)hLJc with ujc = \e\B/{m*c). 
The high degeneracy of the energy levels in absence of a 
potential is associated with a great freedom in the choice 
of a basis of eigenstates for the free Hamiltonian. To diag- 
onalize Hamiltonian ([T]) for an arbitrary potential energy 
landscape V^(r), a very convenient basis^ turns out to be 
the overcomplete set of so-called vortex wave functions 
given by 



«'m,R(r) 



(r|m,R) 

1 { z- Z 



(2) 

. _ |2|2 + |Z|2-2Z2* 

e , (3) 



with z = X + iy and Z = X + iY. The continuous vari- 
able R = {X, Y) constitutes the quantum analog of the 
semiclassical guiding center discussed in the introduction. 
Here we have expressed the wave functions ^ in the 
symmetrical gauge A = B x r/2. Besides being eigen- 
states of the free Hamiltonian, the set of wave functions 
Q has the coherent states character with respect to the 
continuous (degeneracy) quantum number R, which also 
corresponds to a "vortex" -like singularity for r = R. De- 
spite being semiorthogonal, the set of quantum numbers 
|m, R) obeys the completeness relation 



The retarded Green's function (with plus sign in Eq. 
([5])) represents the response of the system to an impulse 
excitation, while the advanced Green's function (with 
minus sign) corresponds to a source wave with a deltalike 
response. Note that the correspondence between Green's 
functions in Eq. ([6]) is nonlocal with respect to the Lan- 
dau level index to, as expected, but quasilocal with re- 
spect to the vortex position R. 

Equation ^ for the electronic Green's function then 
maps^i exactly onto the following Dyson equation for 
the vortex Green's function f/m;m/(R, cj) (from now on, 
we do not specify that the Green's function depends on 
Lu in order not to burden the expressions): 



-co +CO / . X 2fc 

V2 



m"=0 k=0 



X 1 (ax - idy)" Vm;rn"iR) [Ox + idyf ff™'4'(R)-(7) 



The matrix elements Vm;m' (R) of the potential V in the 
vortex representation can be evaluated exactly for an ar- 
bitrary potential provided that the latter is smooth, i.e., 
infinitely differentiable, which is the case for any physi- 
cal potential. They take the form of a series expansion in 
powers of the magnetic length (see Ref. [13 for the detail 
of the derivation) 



(4) 



thus allowing oneSS to use the vortex representation in 
a Green's function formalism, providing unicity of the 
development, related to the analyticity of the disorder 
potentiali^ Note that, however, the nonorthogonality of 
the states prevents to build a perturbation theory solely 
based on wave functions to deal with the potential term 
V{v). 

We have shown in a previous work^i that the elec- 
tronic Green's function associated with Hamiltonian ([T]) 
and satisfying the evolution equation in the energy [uj) 
representation (we set from now on h = 1) 



{lo~H± iS) G^"^(r, r', iv) = S{r - r') 



(5) 



can be written exactly in terms of vortex wave functions 
^'m,R(r) as 



J=0 
3 

E 

fc=0 



V2. 

(to + fc)! 5rn+k,m'+j-k 

Vto! to'! k\{i-k)\ 



(8) 



x{dx+idYf{dx~idYy-^V{K). (9) 



B. Systematic magnetic length expansion 

One method adopted in the paper— in order to solve 
Eq. ([7]) is to search the function g^^ill) under the form 
of a series in powers of the magnetic length Ib , similarly 
to the matrix elements of the potential. 



+ OC' 

5^;™'(R) = E 

3=0 



V2 



(10) 



1^ E *:n'3(r')*™.R(r) 

^ ni.m' 



p=o ■ 



(6) 



Here Ar means the Laplacian operator taken with re- 
spect to the vortex position R, and the term 5 in the left- 
hand side of Eq. ([5]) is an infinitesimal positive quantity 
encoding the boundary condition for the time evolution. 



The functions g'^^-m' (R) ^'^'^ i^&n obtained by solving Eq. 
d?]) order by order in powers of Ib ■ This leads to a closed- 
recursive relational for the functions g'^m-m' (-R) ' which al- 
lows one in principle to obtain an explicit expression for 
9rn-m' oi'dcT j from the knowledge of all other com- 

ponents with subleading order i < j. The leading order 
component can be readily obtained and reads 
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R,AiO), 



R) 



Em~V(Ii)±i6' 



(11) 
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Inserting expression (|lip in Eq. ([6]) and keeping only 
the Ib zeroth order term coming with p = yields the 
compact expression for the electronic Green's function 

X5^';^(°)(R), (12) 

which is a quite simple and general functional of V{R). 
Sublcading corrections up to order were explicitly cal- 
culated in Ref. [2ll . 

C. Limitations of the strict Ib expansion 

Because the term v^^.^, (R) = T^(R) 5m^rn' in the series 
expansion ([8]) is the dominant one for the matrix elements 
of a smooth potential with characteristic length scale ^ ^ 
Ib , one could naively expect that the leading component 
g{o) in (|lip is also the dominant one in the expansion 
(fTO|l for the Green's function. As noted in Ref. [2l|, this 
conclusion has however to be contrasted since it does not 
take into consideration the fact that the (Jb/O expansion 
generates at higher orders systematic terms which may 
be highly singular in energy due to their multiple pole 
structure. This is most clearly seen from Eq. ([6]) for the 
electronic Green's function at coinciding points r = r' 
obtained with the leading order vortex propagator .g'"': 

G^'-4(r,r,c.) (13) 




where integrations by parts and the property 
|*m,R(r)P — |*m,r(R)P wcrc uscd to get the last line 
of Eq. (fT3|l (we have noted above uum — oj ~ Em ± iS). 

Now clearly the truncation of the above Eq. (fT5|) to the 
first p = term is only vindicated provided the integral 
varies on length scales larger than I b , which is not always 
guaranteed, as the vortex wave function spatially extends 
precisely on the scale Ib- If these corrections become 
important, not only the whole sum over p above must be 
kept, but also all terms of similar form that appear within 
the complete vortex Green's function gm;m' (i-e., to all 
orders in Ib)- Let us see what kind of terms one should 
then consider. By inspecting the second line in (jl3[) . one 
is in fact looking for corrections in the vortex Green's 
fmiction at order of the type 

,2. 1 _ 2/||VRyp i%ArV 
^c.,,„-y(R) (c.,„-y(R))3 + (^„,-y(R))2- 

(14) 



Such terms with multiple poles, which indeed start to 
appear in g*^^) (see Ref. [2lj for a complete derivation), 
proliferate at all orders of the Ib expansion, similar to 
the further contributions associated to values of p > 1 
in Eq. ([13]). These corrections to the Green's function 
will not be perturbativcly small whenever one probes 
energies or temperatures smaller than the first charac- 
teristic energy scale appearing above, namely, /s|Vr1/|, 
in which case the leading expression p2)) breaks down. 
Equation ([T^ is however hinting at how a controlled cal- 
culation can be performed: provided that a hierarchy of 
energy scales Zb|VrF| » 1%ArV S> . . . can be estab- 
lished, a systematic resummation to all orders in Ib of 
potential gradient terms will push the validity of the cal- 
culation down to the smaller scale /^ArF, and so on 
and so forth. This idea is quite analogous to the usual 
resummation of classes of Fcynman graphs in standard 
perturbation theory and constitutes the basic motivation 
for the computations that will follow. 



D. Dyson equation in the absence of Landau level 
mixing 

To present the method of resolution of Dyson equation 
(I?!) to infinite order in the I b expansion, we shall focus for 
simplicity on the limit of vanishing Landau level mixing, 
i.e, — > oo with Ib finite. In this case, one can easily 
check that the vortex Green's function becomes purely 
diagonal, 5,„;,„'(R) = gmi'R)Sm,rn' , so that Eq. (O gets 
simplified into 

-l-oo y , \ 2fc -. 

X {dx " idY)" Vm{B.) {dx + idyf 5™^^(R), (15) 

with 

-"W-E^(|ah)V(R). (16) 

Equations (|15p and (|16p are exact in the limit lJc ^ oo 
and valid for any (differentiable) potential V{X,Y). For 
the specific case of a quadratic potential, only the first 
terms k = 0,1,2 and j = 0, 1 of the series appearing, 
respectively, in Eqs. psp and (|16p remain giving rise to 
a nontrivial second-order partial differential equation to 
be solved in Sec. IIIIl Let us first continue considering 
a generic potential and try to simplify at maximum this 
Dyson equation. 

In order to solve Eq. (|15p . it appears very convenient 
to introduce modified vortex Green's function through 
the following change in functions [an insight suggested 
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by the form of the electronic Green's function ([6])] 



-\-oc 



p=0 ^ 



+ 00 



4ar) 



.(R) 



R 



y(R), 



E 

J=0 



(_l)p-j 2J'(TO + j)! 
(p-j)! TO!(j!)2 



(17) 



(18) 



(19) 



where expression (|16p for Um(R) was used to obtain Eq. 
(|19p . After some standard manipulations presented in 
Appendix [Aj one gets the following very compact form 
of Dyson equation (valid for an arbitrary potential, in 
the limit lOc = 00 with Is finite): 

(cc>-S„±z<5)g^'^(R) = 1 



--)vm{R)9^^{R), (20) 



where the notations dx and dy mean that these spa- 
tial derivatives act on the function Vm (R-) only (similarly 
for ^m(R)). Interestingly, and in contrast to the ini- 
tial Dyson Eq. this differential operator starts now 
at order l%dxYvdxY9 (once Dyson equation has been 
properly symmetrized by taking its real part, see Ap- 
pendix so that the change in hmctions (fT7|) and (fT5|) 
manages in principle to perform the whole resummation 
of potential gradient terms to all orders in Ib (this will 
be discussed in more detail in Sec. lIVp . 

Before considering the solution of the transformed 
Dyson equation ((20|) . we need to examine the change 
brought by the mapping ([TT]) in the electronic Green's 
function ([5]), now diagonal in the Landau level index m: 



G^'^(r, r',w) 



+ 00 



^ m=0 



xe 



-^Af 



[*™,R(r')*™,R(r)] , (21) 



where the factors e~^^^ and e'^^^ were combined to- 
gether, and integrations by parts were performed. The 
last step, performed in Appendix [Bl is simply to compute 
the action of the exponential operator in Eq. (PT|) onto 
the product of two vortex wave functions, which finally 
reads: 



where 



[*™,R(r')*m,R(r)] 



.- R[(r 



(R 

— ■)x 



xe 



i*m(R)r 



-A,RVi| 



(22) 



(23) 



with As = (1 - s)/(l -I- s). Form ([22]) will be particu- 
larly useful for subsequent calculations in Sec. IIVI and in 
Appendixes [C] and [El 



III. SOLVING DYSON EQUATION 

A. Absence of curvature: case of an arbitrary ID 
potential or a locally flat disordered 2D potential 

Dyson equation (PU)) . also in its explicit form [Eq. 
(|A7p ]. has the remarkable property that the differential 
operators involve necessarily derivatives of the potential 
in two orthogonal directions. For a ID potential along 
the X direction, the function Vm{X) depends on a single 
coordinate, so that Dyson equation for g^^{X) becomes 
completely trivial and its exact expression (in the limit 
LOc 00) reads 



1 



Vm{X) ± iS' 



(24) 



with Vm{X) defined above in Eq. playing the role of 
an effective potential energy. 

To benchmark expression (j24[) for the modified vor- 
tex Green's function, we consider the exact solution for 
the electronic Green's function that can be derived us- 
ing a standard wave function formalism in the case of 
a parabolic ID potential and prove in Appendix [Cl that 
both approaches lead to identical expressions. This es- 
tablishes that formula (|2ip . with even the lowest order 
vortex Green's function ^m(R), contains the edge states 
physics, which plays an important role in the understand- 
ing of transport properties observed in the quantum Hall 
effect regime<2^i21i^ 

From the present analysis of an arbitrary ID potential, 
one can already guess (see Sec. IIVI for more details) that 
the differential operators appearing in Dyson equation 
((20)) mainly play a role in the case of 2D equipotential 
lines that present a certain amount of curvature at the 
scale of the magnetic length. For a disordered 2D po- 
tential 1^(R), this can occur, e.g., in the vicinity of its 
critical points Rc characterized by Vy(Rc) = 0. For 
an arbitrary smooth potential, and far from its critical 
points, the equipotential lines are locally straight at the 
scale of ^B, so that the modified vortex Green's function 
9i 



''^'"^(R) will be well approximated by the expression 



~9^\R) 



1 



E,^~V{K)±i5' 



(25) 



Once inserted in the electronic Green's function (|2ip . this 
simple result gives the approximate expression 



G«'^(r,r',cj) 



E 



[*™,K(r')*™,R(r)] 



27r?2 

B m=0 



S„ - 1/(R) ± i5 



(26) 



7 



that was proposed with httle detail in our previous Ref. 

Considering that for a smooth 2D potential the equipo- 
tential lines are locally straight on the scale Is (this 
requires a sufficiently large local radius of curvature), 
one can then perform in principle the integration over 
the variable parametrizing distance along the constant 
energy "surface" , T^(R) = const, in the same way as 
done explicitly in Appendix [C] for a pure ID potential. 
One then recovers from the obtained Green's function 
expression the property that the wave functions are lo- 
cally well approximated by translation-invariant Landau 
states with drift velocity c'VVxz/ (|e|i3), as argued in the 
seminal paper by Trugman.— Expression (j26p is however 
quite powerful, because it does not rely on a particular 
parametrization of the equipotential lines, which can be 
cumbersome for a disordered potential, and can be used 
easily by numerically or analytically performing the in- 
tegral over the vortex coordinate R. 

However, as stressed before, approximation (pS)) breaks 
down in the vicinity of the critical points of the potential, 
where the drift velocity locally vanishes. This requires to 
include in the analysis the second-order derivatives of the 
potential V in order to lift the degeneracy of the Landau 
levels, leading to strong quantum effects (quantization 
and/or lifetime), as we will discuss from now on. 



B. Green's functions including curvature effects: 
case of a 2D quadratic potential 



To investigate curvature effects and to determine more 
precisely under which conditions approximation (j25p is 
valid, we expand the arbitrary potential V^(R) around a 
given point Rg, up to quadratic order. This expansion 
appears to be sufficient provided that the gradient and 
the three possible second-order derivatives of the poten- 
tial (locally) never vanish simultaneously, a realistic as- 
sumption. We thus write 

V{R) = V^(Ro) + (R-Ro)- Vr„F(Ro) 

+i[(R-Ro)- VrJV(Ro). (27) 

Inserting expression (|77|) into formula (|18p. we get 



2 ) ^i^^Ir=Ro 



(28) 



«™(R) = V{R) + ^ 



From the symmetrized form ()A7|) of Dyson equation (pO]) . 
we then find that the function g„j(R) is dictated by the 
second-order partial differential equation 

I [{d'yV) dl + {dlV) dl - 2 {dxdyV) dxdY\ ~9l\R) 
-^{uj- E„, - irn (R) ± iS] g^'^{R) = 1. (29) 

The antisymmetrized Dyson equation (jA8|) yields, on the 
other hand, the extra constraint 



indicating that the function g,^j'"^(R) necessarily pos- 
sesses the same equipotential lines as ^(R). We thus 
write 5^'^(R) = /^'^ [EiR)] where £;(R) = V{R) - 
V{Ro) and substitute this expression into Eq. ([^^ to 
obtain a simple ID differential equation obeyed by the 
hmction f,^{^iE): 



E ± tS] /, 



R,A 



C^\E) 
(E) = 1, 



with 



Em — Vm{Ro), 



7 



dj^Vd^V 



{OxdyVy 



2 



■0 



'2dxVdYVdxdYV 



R, — R,n 



(31) 

(32) 
(33) 

(34) 



The coefficient 7 is nothing but the determinant of the 
Hessian matrix of the potential V, with a prefactor l%/4:. 
Its sign determines the nature of the critical points Rc at 
which IVr^I vanishes. A saddle point is characterized 
by 7(Rc) < 0, while a strictly positive 7(Rc) indicates 
the presence of a local maximum or minimum. 

Differential equation (|3ip can be solved by Fourier 
transforming to time, see Appendix [D] so that the 
Green's function is given by the explicit formula 



s^'^(R) 



with 
^m^(Ro,i) 
r{t) 



dth^^{Ro,t)e 



-i[y(R)-y(Ro)]r(t) 



(35) 



cos(/7i) 

— tan {^t) . 
\/7 



-i(l]/'()T{t)+i{uirn+V 



(36) 
(37) 



{dxV)dY~g'rf,'''-{dYV)dx~g, 



R,A 







(30) 



Noticeably, when 7 > the function r(t) is a periodic 
function of time t, so that Green's function gm must dis- 
play discrete poles, and quantization of energy levels in a 
confined potential is recovered (this is further discussed 
in Appendix IE|). We stress that such success of the vortex 
formalism was far from granted, because one has started 
with a family of wave functions labeled by the continuous 
quantum number R. 

For 7 < 0, the functions cos and tan in Eqs. ([55]) 
and (j37[) are to be replaced by their hyperbolic coun- 
terparts cosh and tanh, respectively, so that the kernel 
1/ cosh(-y— 7t) obviously introduces lifetime effects in the 
description [the convergence of integral ([55)) over the time 
is now ensured by this term and no more by the cutoff 
function exp{^5t), as was the case for 7 > 0]. Clearly, 
the vortex self-energy obtained from Eq. (j35p displays 
an elastic scattering rate proportional to ^—7, a clear 
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signature of quantum tunneling at saddle point with im- 
portant consequences for transport properties (see dis- 
cussion in Sec. |V|. This allows us to make the crucial 
physical identification between scattering mechanism and 
negative curvature of the potential in the quantum Hall 
regime. 

It is interesting to note that the strong quantum effects 
(quantization or lifetime) exhibited by the exact quantum 
solution ([35]) are dictated by the quantitity .^7, which in- 
volves the square root of the second-order derivatives of 
the potential. Clearly, they thus can not be fully cap- 
tured via a finite expansion in powers of the magnetic 
length which can only generate integral powers of the 
derivatives of the potential; see Sees. Ill Bl and III CI This 
impossibility to approximate quantum effects at finite I b 
in a controllable way with the I b expansion illustrates its 
asymptotic character. 

The function /i^''^(Ro,t) depends on the reference 
point Ro via the coefficients 77 and ojm for a generic 
quadratic potential, and possibly also via the coefRcicnt 
7 for a potential characterized by higher derivatives. The 
geometric parameters 7 and rj are basically small coeffi- 
cients for a potential ViTl) which varies smoothly at the 
scale Ib- If we literally take 7 = ry = 0, we find again 
expression ([25|) for the function gm(R). We will show fur- 
ther in which circumstance it is nevertheless required to 
keep the dependence on the coefficients 7 and/or jy in the 
Green's function to correctly describe the local physical 
observables. 

Making use of expression ([55]) together with Eq. pT|) . 
we obtain that the electronic Green's function reads fi- 
nally 



G-«'^(r,r',w) 



d^R if? 



E 



[«'™,R(r')*™,R(r)] 



jdt h^^^iH^, t) ^-^[Vm-V(Tlo)]r(t)^ (38) 



Expression (|38p is the main mathematical result of this 
work. It is exact in the limit Wc — > 00 for any quadratic 
potentials. In particular, it holds for quadratic confining 
potentials simulating closed systems, such as quantum 
dots, as well as for nonconfining quadratic potentials cor- 
responding to open systems, such as quantum point con- 
tacts. Physical implications of this result are discussed 
in Sec.|Vl while further mathematical simplifications will 
now be performed in order to extract relevant physical 
observables. 



IV. LOCAL DENSITY AND CURVATURE 
EFFECTS 

A. Simplifying the Green's function expression 



Expression (|38p for the Green's functions can be calcu- 
lated further in different ways. One possibility is to use a 
paramctrization of the equipotential lines of ViTC). Such 



an approach appears, however, not very practical for a 
generic random potential. Actually, it turns out that 
the two-dimensional integral over the position R can be 
performed analytically when V^(R) is expanded up to its 
second derivatives around the point Rq. For a quadratic 
potential, the Green's function can be rewritten at the 
final stage as a single one-dimensional integral over the 
time variable i, as will be shown in this section. For the 
numerics, this appears to be more easily tractable than 
a direct computation of formula (|38p . 

Note that for a quadratic potential, formula ([55]) is ac- 
tually independent of the choice of Rq. This can be eas- 
ily checked by taking the gradient of expression p8|) with 
respect to Rq and considering that, besides the explicit 
term y(Ro), the dependence on Rq is also contained in 
the function /im(Ro, t) through the coefficients ?/ and ujm 
(the other coefficient 7 is independent of Rq in the par- 
ticular case of a quadratic potential) . The independence 
of the electronic Green's function then follows from the 
relation Vrq?? = 7 Vr(,V(Ro). 

For a smooth arbitrary potential 1^(R), result ([55)1 is 
expected to give a very good approximation to the elec- 
tronic Green's function provided that the temperature 
exceeds the energy scales associated with the third order 
(and beyond) derivatives of the potential. Contrary to 
the case of a quadratic potential, formula ((38)l will now 
depend on the reference point Rq, which thus has to be 
chosen appropriately. The natural choice appears to be 
Ro = (r + r')/2. 

Inserting formulas ([^ and (^5]) into Eq. using 
the expansion (P?]) of the potential F(R) up to quadratic 
order with Rq = (r -I- r')/2 = c and evaluating the re- 
sulting Gaussian integrals over the variable R, we get 



+00 „ 

G«'^(r,r',a;)=5] dt 



exp 



2(4) A3i||W(r,r',t)|^+4j))(r,r',t)r(t) 
"4 A^^+iA^C-r<,t)—yT^{t) 



(1 + s)^Al+iAsCr{t)-jr^t) 



with 



W(r,r',t) = VcV^(c) + z X 



llrit)- 



(39) 

/ s=0 

(40) 



C 



[{W X i) ■ HvlJW X z)] , (41) 

(42) 



|AcV(c), 



where Hv\^ is the 2x2 Hessian matrix composed of the 
second derivatives of the potential V taken at position c 
(the elements of the matrix are given by [Hyjij = didjV). 
For r = r', we have the simplifications W(r, r, i) = 
'VrV{r) and ?)(r,r,t) = 7]{r) [function defined in Eq. 
([M)) ]. Note that for a potential V characterized by 
derivatives of order higher than 2, formula ([55)) yields 
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only an approximate result. In this case, all the geomet- 
ric coefficients, including 7 and ^, depend on the center 
of mass position c. 



B. Local density of states 

We now aim at computing the local electronic density 
defined by 



n(r) = /do; ni?(tj)p(r, cj), 



(43) 



where the local density of states p{r,uj) is directly ob- 
tained from the retarded Green's function at coincident 
positions as 



p(r,u;) = -ilm[G^(r,r,c.) 
n 



(44) 



Here np{uj) = [1 + cxp([ijj — ^]/T)]^^ is the Fermi- 
Dirac distribution function, T the temperature, and /i 
the chemical potential. We thus need the simpler form 
of expression 



G«'^(r,r,c.)= 2 Ai''"^^""'*^ ^ 



27r/| 



exp 



'■{t) AJ%\^,V{v)\'+iir^(v)r{t) 
4 Ai+iA^CT{t)~-tT-'{t) 



as" 



(1 + s)^A^,+lAsCT{t)-^T^{t) 



.(45) 



To simplify further the expression of the local density, 
it is then required to consider the explicit expression (j36p 
for the function h^{r,t) and insert it into Eq. (|45ll . In 
order to do the integral over lo in expression (|43p , we first 
introduce the change in variable lo' = lo ~ fi and decom- 
pose the exponential factor in the numerator depending 
on Lo' as exp(ict''i) = cos{Lu't) + ism{Lo't). The integral 
over the energy to' in Eq. (I43p coming with the first term 
cos^Lo't) is then performed by writing the Fcrmi-Dirac 
distribution function as 



nF{uj' + fj.) 



1 — tanh 



2T 



(46) 



On the other hand, the second contribution to the inte- 
gral coming with the term sm{Lo't) is calculated by 
using the result 



, sin(w'i) ttT 

OLO 



1 + e"'/T 



sinh(7rTi) ' 



(47) 



Finally, we find that the local density takes the form of 
a simple integral over the time t 



s=0 



n(r) 



27r/| 



T 



- + Im V / 

2 ,^0-^0 sinh[7rTt] 



,^t[^L-Er„-(m+l/2■K-V{r)]^^lt~r{t}] 



1 9™ 
m! ds^ 



exp 



cos(77<) 

4 Ai+iA,CT{t)--iT^(t) 



(1 + s)^A^,+zA,CTit)~-fr^t) 



(48) 



This formula is exact for any quadratic potential in the 
absence of Landau level mixing. To illustrate this strong 
statement, we prove in Appendix [E] its equivalence with 
the expression for the local density that can be derived 
by standard means in the specific case of a circular 2D 
parabolic confinement (note that we have already shown 
the correspondence in the different case of a ID parabolic 
potential at the level of the Green's functions, see Ap- 
pendix [C|. This shows that quantization effects, i.e., the 
presence of a discrete energy spectrum, are fully captured 
in the vortex representation, despite not being completely 
explicit in formula (|48p . The latter equation has thus a 
relatively general character since it contains under a com- 
pact and unified form the cases of confining and noncon- 
fining quadratic potentials. Note that expression (|48)) is 
naively problematic for the saddle-point quadratic po- 



tential model because the energy spectrum in this case 
is unbounded from below, but relative density variations 
arc, on the other hand, perfectly well defined. 

Of particular interest is the derivative of the local den- 
sity with respect to the chemical potential which can be 
directly probed by the differential tunneling conductance 
in a scanning tunneling spectroscopy (STS) experiment 
(provided that the tip density of states is constant in the 
studied energy range): 

p«^«(r,/i,r) = ^ = j du[-n'^{uo)]p{v,u). (49) 

At zero temperature, this yields the local density of 
states at the chemical potential energy, p(r, /j,), since then 
—n'p{uj) ~ 5{uj — [I). Using formula we directly get 
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Tt 



+ 00 (. + 00 

(it 

-JO sinh[KTt\ 



gjt[M--E™-(m+l/2)C-V(r)]gi2M[t-r(t)] 



1 9'= 



cos(^t) 
exp 



T^jt) Ajl\-V,V{r)\^+4iri{T)T{t) 
4 Ai+iA^C'r{t)-jT^{t) 



■m!9s™ I (l + s)^A2+a,Cr(i)-7^'W 



(50) 



Contrary to the local density formula, expression (j50p 
is well defined for the saddle-point quadratic potential 
model because it involves only states around the energy 
/i. Formula (|50|) for the local density of states is exact 
for any quadratic potential. One may wonder about its 
accuracy for an arbitrary potential landscape which is 
smooth on the scale of the magnetic length. We shall in- 
vestigate this question by a careful quantitative analysis 
in the next subsection. 



C. Quantitative aspects: when do gradient and 
curvature corrections need to be included? 



In order to illustrate on a concrete example how suc- 
cessive steps in the resummation of leading derivatives 
of the potential really operate, we focus here on the 2D 
circular confining potential 



Vir) 



1 



(51) 



whose explicit solution is given by the so-called Fock- 
Darwin states (see Appendix |E]), and investigate the 
temperature-dependent local density of states 

The simplest approximation scheme, which amounts 
to view the potential term (|51|) in a purely local manner, 
i.e., V{R) ~ V{Ro), is obtained by setting |VV^| = C = 
7 = = in Eq. ([50)1 . This obviously recovers the usual 
semiclassical guiding center result: 



00 p+oa 



n"'0 



'dt Tt cos(t[^ -Em- V{r)]) 



m=0 " ^ 



sinh [irTt] 



2ttI 



2 +00 

-^Y.n'AEm + Vir)]. 



(52) 



m=0 



This result is in fact accurate as long as one considers 
temperatures higher than the energy scale associated to 
the drift motion, namely, ^s|VRy|. At lower tempera- 
tures, the resummation of all leading gradient contribu- 
tions needs to be performed, which corresponds to con- 
sidering the potential as locally flat (in the geometrical 
sense): V{R) ~ V(Ro) + (R- R-o) ■ VRoy(Ro). This 
calculation can in fact be achieved with the previously 



obtained results, setting C = 7 = ''? = 0in Eq. (|50p : 

'-00 ^+00^^ Ttcos{t[^i- Em-V{r)]) 



m=0 " " B 

Qm j exp 



sinh [nTt 

e il;|V.V(r)p 



m! 9s™ 



{l + s)A, 



(53) 



s=0 



Clearly the scale /s| VrI^I provides a cutoff in the above 
integral, so that the single pole divergence associated to 
the derivative of the Fermi-Dirac distribution function in 
Eq. ([52|) is regularized. 
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FIG. 1: (Color online) STS density of states at peak position 
Tpcak, in units of {2-kI%)^^ , as a function of T /TujJc, for circu- 
lar potential (|5ip with uq/uc — 0.1, and ^/huic = 0.8 (lowest 
Landau level). Semiclassical result (|52|l . gradient resumma- 
tion (|53p . and full quadratic solution (|50p including curvature 
effects are presented. 

Fig. [Ddisplays the STS local density of states as a func- 
tion of temperature for a fixed chemical potential and the 
particular position rpcak given by 11 - Em - V(rpeak) = 
0, according to semiclassical expression (|52p. the lead- 
ing gradient approximation (|53p and the exact solution 
([50)1 which includes also curvature effects from the full 
quadratic dependence of potential pT)) . For the sake of 
simplicity, we have considered ^ = O.Sfiuc, which corre- 
sponds to filling the lowest Landau level to = only. 

Clearly in Fig. [1] the semiclassical approximation is 
only valid in a high temperature regime and breaks down 
below the energy scale ^bIVrI/I. The departure of semi- 
classical expression (|5^ from the exact one [Eq. (|50p ] 
is also easily seen in the Figs. [2] to SI where the spatial 
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FIG. 2: (Color online) Similar data as in Fig. [T] now as 
a function of radial distance r/Zs, for the high temperature 
T/fttJc = 0.1. 
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FIG. 3: (Color online) Similar data as in Fig. [T] now as a 
function of radial distance r/ls, for the intermediate temper- 
ature T/hiVc = 0.01. 



saturation at intermediate temperatures [where semiclas- 
sical result ([5^ is already quite inaccurate] of the STS 
density of states at peak value; see also Fig. [3l Decreas- 
ing further the temperature, curvature effects associated 
to the small geometric energy scale ^ begin to be felt; 
see Figs. [T] and IH Departure from the leading gradient 
result is manifest by a final divergence of the exact den- 
sity of states at the peak position in the limit T — s- 0, 
since the Fock-Darwin energy spectrum is discrete (with 
level spacing ,^7) due to the confinement. 

As a final remark, the above discussion is quite instruc- 
tive, as it clearly shows under which conditions curvature 
effects associated with second-order derivatives of the po- 
tential can be neglected, namely, when temperature is 
higher than the energy associated with curvature. There- 
fore, successive approximation schemes can be devised for 
a smooth arbitrary (disordered or not) potential leading 
to controlled expressions for the local density of states. 
The whole scheme is indeed based on the existence of a 
hierarchy of local energy scales of the type l^dl!'V{r) . Ex- 
pression ([50]) . which includes all second-order derivatives 
of the potential, thus provides an accurate estimate as 
long as the temperature is larger than cubic and higher 
order derivatives of the potential. In particular, it is also 
valid near saddle points of the potential landscape, where 
the drift velocity vanishes. It is thus extremely useful for 
interpreting local STS experiments such as in Ref. [2^ and 
completely bypasses the need to diagonalize numerically 
a complicated random Schrodingcr equation. 



V. DISCUSSION: ON THE FUNDAMENTAL 
IMPORTANCE OF OVERCOMPLETENESS 



dependence of the STS density of states is plotted for 
three different temperatures. Expression (|53p which has 
a greater domain of validity than the semiclassical one 
turns out to match the exact result down to temperatures 
of the order of the curvature energy scale The expo- 
nential cutoff in Eq. ([53|) manifests itself on Fig. [T] by a 
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FIG. 4: (Color online) Similar data as in Fig. [T] now as 
a function of radial distance r/ls, for the low temperature 
T/huJc = 0.001. 



In the light of the technical results derived in the previ- 
ous sections, we formulate here some general conclusions 
on very fundamental issues in quantum mechanics such 
as the emergence of classicality and the microscopic ori- 
gin of time irreversibility. 



A. Emergence of classicality in quantum mechanics 

It is well-known that the classical Hamilton- Jacobi 
equations of motion can be derived from the quantum- 
mechanical Schrodinger equation when terms having h 
as prefactor can be disregarded. In other terms, classi- 
cal mechanics is clearly a limit of quantum mechanics. 
However, capturing the precise mechanism responsible 
for the emergence of the classical behavior in the physi- 
cal properties of the system within a fully quantum me- 
chanical framework, i.e., at h finite, appears much more 
complicated. The essential reason is that establishing 
the quantum-classical correspondence requires not only 
to consider the equations of motion but also the states 
of the system. And the limit h s- appears to be much 
more singular for the wave functions than for the energy 
spectrum. When dealing with this limit, we are imme- 
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diately confronted with a conceptual problem relying on 
the fact that quantum mechanics is originally formulated 
in a Hilbert space spanned by a countable basis of square 
integrable states, while classical dynamics occurs in a 
continuous phase space. We are therefore in a delicate 
position to reproduce the basic structure of the classical 
phase space. 

In the particular problem under study in this paper, 
the set of vortex states |m, R) introduces from the very 
beginning in the quantum description a continuous repre- 
sentation for the quantum numbers. Because they obey 
in part the coherent states algebra^^ (note that the vortex 
states are very peculiar coherent states in so far as they 
present the coherent character only via the degeneracy 
quantum number R and not via the eigenvalue quan- 
tum number m, so that they can be also eigenstates of 
the kinetic part of the Hamiltonian, in contrast to fully 
coherent states), and especially a completeness relation 
[Eq. ([4])], we can legitimately use the vortex represen- 
tation for the spectral decomposition of Hamiltonian ([1]) 
provided that the potential F is a smooth functioni^ As 
an original motivation to work preferentially with these 
stateS)2£ the quantum numbers |m, R) provide a very in- 
tuitive and clear physical connection to the classical dy- 
namics for the free Hamiltonian when considering the 
de Broglie-Madelung hydrodynamic pictureS^i^ of the 
Schrodinger equation: the quantization of the kinetic en- 
ergy into Landau levels stems only from the interference 
of the electronic wave function with itself due to the com- 
pletion of a circular orbit around the position R, where 
the phase of the wave function is ill-defined. The price to 
pay for the continuous aspect, i.e., for introducing over- 
completeness into the quantum-mechanical formalism, is 
the nonorthogonality of the states with respect to the 
degeneracy quantum number R, which reflects the quan- 
tum indeterminacy in the positions of the vortices and is 
accounted for in formula (jH) by associating the elemen- 
tary area 2ttI^ to the incremental area in the integration 
over the vortex positions. 

Being better armed to capture the transition from 
the quantum to classical, it is not completely a surprise 
that we find that the vortex representation leads at the 
mathematical level to a systematic and straightforward 
expansionSi in powers of the magnetic length (which, we 
remind, plays the role of an effective Planck's constant 
in the present problem) of the vortex Green's functions, 
and thus of the physical observables. Therefore, it turns 
out that overcompleteness is clearly not a drawback but 
an advantage at the technical level! However, behind 
this mathematical aspect, we also see a very fundamen- 
tal physical aspect, which is rarely considered in quan- 
tum mechanics when choosing a peculiar representation 
of states. Obviously, the vortex representation offers the 
unique opportunity to derive quantum expressions with- 
out having to implement the complete explicit form of the 
potential V. This is exemplified by the exact compact 
formula ([551) f^i' the Green's function which embraces all 
possible cases of quadratic potentials. The generic form 



of this result actually encodes the stability of the vortex 
states. Indeed, the Fock-Darwin states (jE3p which corre- 
spond to the exact eigenstates of Hamiltonian ([1]) in the 
presence of a circular parabolic confinement and have 
the rotational symmetry (see Appendix [E]) appear to be 
very unstable: one can not expect the confinement to be 
perfectly circular under realistic conditions, so that the 
real physical state certainly does not obey the rotational 
symmetry. In contrast, the vortex states which enclose 
no preferred symmetry turn out to be stable with respect 
to an arbitrary small asymmetrical smooth perturbation 
of the potential landscape. From this robustness prop- 
erty, we can expect them to be the real physical states, 
i.e.; the most predictable ones in an experiment. 

Interestingly, we have an illustration with the present 
study for the process of superselection of states put for- 
ward by Zurek^ to explain the emergence of the clas- 
sical behavior from a quantum substrate. The only im- 
portant difference is that we are somewhat accounting 
here for an intrinsic mechanism of classicality. Indeed, 
it is customary in quantum mechanics to appeal to ex- 
trinsic degrees of freedom brought by an environment 
(surrounding the studied quantum system) to explain 
the appearance of classical properties through decoher- 
ence processes. As developed by several authors (see the 
review^!), the environment prevents certain quantum su- 
perposition of states from being observed as a result of 
their high instability. Only states that survive this pro- 
cess of coupling to the environmental degrees of freedom 
have predictable consequences. As shown by Zurek et 
al^ in a model of weakly damped harmonic oscillator, 
coherent states, which are known to be the closest states 
from the classical limit, are minimally affected by the 
coupling to the environment. Due to this robustness, 
they emerge as a preferred set of states. 

In the present problem of the electron dynamics in 
a high magnetic field, we clearly see under which con- 
ditions the overcomplete vortex representation becomes 
effectively selected by the dynamics. Indeed, we have 
noted that formula (|38|) derived in the vortex representa- 
tion reproduces the exact Green's functions in the simple 
integrable case of a 2D circular confining potential (Ap- 
pendix |E| . The system actually does not exhibit yet a 
preference for the overcomplete set of vortex states over 
the complete set of Fock-Darwin eigenstates. In contrast, 
the case of a quadratic saddle-point potential which sim- 
ulates an open system and introduces a dynamical in- 
stability seems quite instructive. Indeed, the conven- 
tional approach of quantum mechanics with square in- 
tegrable wave functions turns out to be inadequate to 
determine the energy spectrum, so that one usually has 
to resort to another formalism, namely, the scattering 
states quantum formalismj2^ These difficulties are man- 
ifestations of the fact that the spectral problem for un- 
stable unconfined dynamical systems is not computable 
in the Hilbert space. The overcompleteness of the vortex 
representation in this specific case of saddle-point poten- 
tial shows precisely its relevance by allowing one to solve 
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the dynamical equations exactly on the same footing as 
in the confining cases in a Green's function formalism. 
One can thus argue that the overcomplete set of vortex 
states is naturally favored by the instability of the dy- 
namics. Noticeably, the basic dynamical object in the 
vortex representation appears to be no more the wave 
function but the Green's function. By inspecting the 
form of the generic Green's function ([55]) . one notices 
that the latter can not be written explicitly as a product 
of two wave functions (as is usually the case when us- 
ing a complete representation) due to the presence of the 
nonlocal operator exp [— (Z^/4)Ar] acting on the vortex 
wave functions [see also Eqs. (|22p - ((23)) ]. This reflects the 
overcompleteness of the coherent states basis with the 
two-dimensional continuous quantum numbers R associ- 
ated to the vortex position. It is therefore clear that it 
is not possible to get an single expression encompassing 
all possible cases of confining and unconfining quadratic 
potentials in terms of wave functions eigensolutions of 
the Schrodinger's equation. This general result can only 
be achieved through the introduction of an overcomplete 
basis of physical states. 



B. Time irreversibility 

An attractive feature is the close links existing between 
the transition from quantum to classical (as a result of de- 
coherence) and time irreversibility. By time irreversibil- 
ity we mean the time asymmetry due to a preferred di- 
rection of time, as shown by decaying states. While 
quantum mechanics is able to provide a clear and suc- 
cessful dynamical foundation to the idea of quantum lev- 
els, the problem of decaying states with lifetimes remains 
somewhat obscure and controversial. These difficulties in 
identifying the physical roots of irreversibility rely essen- 
tially on the fact that the microscopic dynamical equa- 
tions are time reversible, whereas complex macroscopic 
systems are always characterized by a time-asymmctric 
evolution. Consequently, it is generally believed that irre- 
versibility arises from the macroscopically large number 
of degrees of freedom affecting the time evolution of a 
nonisolated systemi^ 

There have been many different approaches to de- 
rive an irreversible dynamical evolution starting from the 
Schrodinger equation. The most popular one^ is to con- 
sider the microscopic (integrable) system as a part of 
a larger Hamiltonian system which has many degrees of 
freedom (the environment or reservoir). Then, after trac- 
ing over the environmental degrees of freedom (the latter 
are disregarded because uncontrolled and unobserved), 
the dynamics of the (open) quantum system is no more 
described by the Schrodinger equation, which is expected 
to be applicable only to a closed system. Other possibil- 
ities are to solve quantum-mechanical equations by deal- 
ing directly with tractable models of the environment, 
such as the consideration of a collection of harmonic oscil- 
lators. The common denominator of all these approaches 



is to associate time asymmetry with the external influ- 
ence of a reservoir or a measurement apparatus. Irre- 
versibility thus seemingly has an extrinsic root. 

In order to better clarify its possible link with the in- 
herent dynamics of the system, Prigogine et aZi ^^'^^i'^^ de- 
manded that irreversibility be rather directly connected 
with the Hamiltonian of the microscopic quantum sys- 
tem, in spite of introducing extra dynamical assumptions 
(because, after all, the division of a global system into a 
system and an environment is artificial and rather a mat- 
ter of taste). These authors^^i^ii^ used extensions of the 
traditional Hilbert space through the introduction of a 
nonunitary change of representation and argued with a 
few simple examples that time asymmetry may sponta- 
neously arise in systems whose dynamics is nonintegrable 
in the Hilbert space of quantum mechanics. Then, the 
problem of integration and irreversibility both enjoy a 
common solution in the extended space. 

In the Hilbert space quantum mechanics, the time 
evolution described by the Hamiltonian must be time 
reversible, leading to a widespread belief that intrin- 
sic irreversibility simply docs not exist. Moreover, for 
the nontrivial physically interesting systems, the com- 
putability of the spectral problem is generally limited, the 
state of the art offering only perturbative and/or effec- 
tive approximate solutions. In such systems, irreversibil- 
ity does appear in the derivation, but as the result of 
supplementary approximations to the Hamiltonian for- 
malism of quantum mechanics. A well-known example 
in condensed-matter physics is the case of a disordered 
system for which elastic lifetimes in the spectrum are 
obtained by averaging over disorder configurations^^ In 
brief, in order to clarify an intrinsic mechanism of ir- 
reversibility, it is of valuable interest to find nontrivial 
physical systems which are sufficiently simple to allow 
exact time-asymmctric solutions. 

We strongly believe that the exact solution for the elec- 
tron dynamics in a high magnetic field and a given yet 
arbitrary quadratic potential presented in this paper pre- 
cisely offers such an opportunity. We have noted in Sec. 
mil that the Green's functions are characterized by the 
presence of lifetimes in the case of saddle-point poten- 
tials (when the geometric curvature 7 < 0), meaning 
that time symmetry is broken. We thus obtained irre- 
versibility without appealing to extra dynamical consid- 
erations, such as an environmental coupling. In other 
terms, we are basically in the scenario depicted by Pri- 
gogine et aljS^'^^j^ 

One may naturally wonder how the time-reversible 
Schrodinger equation can then lead to irreversible pro- 
cesses at the mathematical level. It is often believed that 
the complex poles of the Green's functions correspond 
to eigenvalues of a non-Hermitian operator. In contrast, 
we would like to point out that a broken time symme- 
try exhibited by the states is not necessarily in contra- 
diction with a time-invariant Hamiltonian if a mathe- 
matical theory is used that makes a distinction between 
states and the Hermitian Hamiltonian operator. Actu- 
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ally, the dynamics remains here time symmetric but is 
realized through an overcomplete representation which 
permits a broken time symmetry for the states. A com- 
plete (countable) representation for its part does not al- 
low time-asymmetric solutions. The overcomplete vor- 
tex representation provides a more general type of spec- 
tral decomposition of the Hamiltonian operator, which 
is merely based on the use of Dirac's bra-kct formalism. 
The extension of the eigenvalue problem to the complex 
plane is then purely a qualifying feature of the instability 
of the dynamics, thus revealing an intrinsic irreversible 
character of the evolution of the states. 

It has been stressed by several authors^^^i^i^ that 
the natural setting of quantum mechanics is the rigged 
Hilbcrt space rather than the Hilbert space alone. The 
rigged Hilbert space is just an extended space consisting 
of the Hilbert space equipped with distribution theory 
and was originally introduced into quantum mechanics 
to give a mathematical justification of Dirac's bra-ket 
formalism. It establishes rigorously that the spectral de- 
composition formula acquires meaning in the continuous 
spectrum as well as in the discrete spectrum, and al- 
lows the appearance of complex eigenvalues. Plane-wave 
eigenvectors, which are generalized eigenvectors in the 
space of tempered distributions, are basic examples of 
these elements of the rigged Hilbert space which do not 
live in the Hilbert space. They are routinely used in the 
scattering states formalism, which contains an arrow of 
time hidden in the choice of time asymmetric boundary 
conditions: The consideration of in- and out- plane wave 
states asymptotically far from the scattering region is in- 
deed a statement of causality expressing the fact that the 
state at a given position is determined by the action of 
a source at a retarded time. Note that causality is natu- 
rally accounted for in the definition itself of the retarded 
and advanced Green's functions. However, in this case, 
the presence of the infinitesimal quantity S in the dynam- 
ical equations [see Eq. ([5])] does not automatically imply 
a broken time symmetry for the physical states. For this, 
one needs in addition to have a dynamical instability oc- 
curring in an unconfined system, i.e. scattering events. 

The introduction of a continuous ingredient plays an 
important role in all microscopic derivations of irre- 
versible processes. With the consideration of asymptotic 
in- and out-plane wave states, the scattering formalism 
presupposes the existence of a continuum via the pres- 
ence of reservoirs, so that irreversibility finally acquires 
within this approach an extrinsic character. Moreover, 
this formalism is specifically limited to open systems, 
thus antagonistic to the Hilbcrt space quantum mechan- 
ics of closed systems. In this paper, we have shown 
that, by using an overcomplete representation of coherent 
states belonging to the Hilbert space such as the vortex 
states, it is possible to embed quantum theory in a wider 
formalism of which Hilbert space quantum mechanics of 
closed systems would become a special case. Moreover, 
in this approach quantization effects and lifetime effects 
are naturally treated on the same footing. The contin- 



uous ingredient is contained into the overcompleteness 
property of the chosen set of quantum numbers. As a 
price to pay when working in a coherent states represen- 
tation, it requires giving up the wave functions as the 
fundamental quantity of quantum theory and replacing 
them by Green's functions. It is worth emphasizing that 
the overcompleteness does not necessarily imply a loss 
of information and time symmetry breaking. For this, 
we need in addition an instability of dynamical motion 
related, e.g., to the presence of saddle points in the po- 
tential landscape. In this case, the overcompleteness of 
the representation^ is necessary to obtain a solution of 
the spectral problem. The basic reason is that the cross- 
ing of the equipotential lines at the saddle-point energy 
(which schematically looks like a collision process and can 
be seen as a bifurcation of a path) together with the ope- 
ncss of the system destroys the trajectory as well as the 
Hilbert space description. Therefore, the phenomenon of 
instability somehow imposes to deal directly with proba- 
bilities to describe the dynamical evolution of the physi- 
cal states (which necessarily belong to the Hilbert space). 
It is worth noting that we then obtain a description which 
from the point of view of its structure is isomorphic to 
classical mechanics. 

We have seen that irreversibility arises as a selec- 
tion principle from the time-invariant Hamiltonian. The 
states selected by the unstable dynamics appear thus 
to be less symmetric than they would seem to follow 
from the Hamiltonian description. This situation is ac- 
tually reminiscent of the well-known spontaneous sym- 
metry breaking as it occurs in ferromagnetism. In the 
presence of a dynamical instability, bra and ket vortex 
states describe just physically distinct states. Finally, we 
note that a critical ingredient to obtain the time symme- 
try breaking in our solution is to consider quantum tun- 
neling within an infinite system, i.e., unconfined spatially 
(otherwise, the physical quantum numbers describing the 
dynamics are necessarily discrete and the evolution uni- 
tary). 



VI. CONCLUSION 

In this paper, we have built a Green's function formal- 
ism based on the use of an overcomplete semicoherent 
vortex representation to study the electron quantum dy- 
namics in high magnetic fields and in a smooth potential 
landscape. Within this formalism, we have shown that 
it is possible to derive in a controllable way approximate 
quantum expressions, e.g., for the local density of states, 
for an arbitrary potential smooth at the scale of the mag- 
netic length. Moreover, we have obtained in the limit of 
negligible Landau level mixing an exact expression for the 
electronic Green's function which encompasses all possi- 
ble cases of quadratic potentials. We have argued that 
this generic result, which is rendered possible by the use 
of an overcomplete representation of states belonging to 
the Hilbcrt space, is a manifestation of a stability prop- 



15 



erty of the vortex quantum numbers. We have shown 
that the overcompleteness feature of the vortex repre- 
sentation does not introduce de facto a loss of informa- 
tion, since we are able to reproduce the solutions for the 
exactly solvable (integrable) cases of parabolic ID and 
2D confining potentials, which can be obtained by stan- 
dard wave function calculations. In contrast, wc have 
found that a loss of information, associated with the in- 
troduction of a probabilistic description of the physical 
processes, and concomitant with the appearance of life- 
times (synonymous of time symmetry breaking) , arises in 
the saddle-point quadratic potential model. The vortex 
representation turns out to be especially relevant in this 
latter case of quadratic potential by providing in the limit 
of negligible Landau level mixing exact physical insight 
into the quantum tunneling processes originating at the 
saddle point. Therefore, we have explicitly proved that 
time irreversibility does not result from supplementary 
approximations to the Hamiltonian formalism of quan- 
tum mechanics, but just naturally arises in the spectral 
decomposition of the Hamiltonian from the formulation 
of dynamics in this overcomplete vortex representation 
of states. With the present analysis, we deduce that 
the minimal necessary ingredient to get solutions from 
the Hamiltonian formalism which exhibit a broken time 
symmetry is to have an instability of the single-particle 
dynamics, as occurring from quantum tunneling at the 
saddle points of the potential landscape, which manifests 
itself in an unconfincd (thus open) system. Therefore, 
besides permitting to capture the transition from quan- 
tum to classical in an efficient way, the overcompleteness 
property of the representation allows the introduction of 
an intrinsic irreversibility on the microscopic level. 



APPENDIX A: DETAILS ON THE MAPPING OF 
DYSON EQUATION IN THE HIGH FIELD LIMIT 



Dyson equation (jl5[) has been rewritten in the ujc = oo 
limit and we aim here at getting a simpler yet equivalent 
form that trivializes the problem of local potential gra- 
dients. This can be achieved through the substitution of 
functions (flT)) and (fTS]) . which clearly gives 



(c^-i?„±i<5)5^'^(R) = l + 

2k 



(Al) 



Going to Fourier space permits to rewrite the right-hand 
side of expression (|A1[) as a single global operator. In- 
deed, defining 



gm{R) = jd\~gr,M)e'^-^, (A2) 
^;„(R) = Jd^piMe'P-^, (A3) 



Acknowledgement 



T.C. acknowledges interesting discussions with D. M. 
Basko. 



and inserting these expressions into the right-hand side 
of Eq. (jAl[) . important simplifications occur: 



+00 / . s 2fc -, 
i-B \ 1 



^'q^'pE [^) ^vM-gM [{^P.+Py)i^q. - qyt eMiP+^f-p'-^^^iP+^y^ = 

d2qd2p{}^(p)g^^(q)eT^(''p-+P«)(''9--'J«)e%-P-'ie*(P+'i)-^ = 

d2qd2pi)™(p).g™(q) e''4-(Py9--Px9y)e*(p+q)-R-. (A4) 

I 

The global operator e^^BiPvi^-p^iv)/"^ above can then be We note in passing that the other Dyson equation ( i.e., 
written back into real space, providing the final expres- G = Gq + GVGq) provides a second equation satisfied by 
sion given in Eq. 
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oo / , \ 2k 



the function gm'- 

X 1 (9x + idyf «m(R) {dx - idvf ff^'^(R), (A5) 

which may be mapped in a similar way onto the following 
equation for the function 

{L0-E^±i5)~g^^{Ii) = l 

A more explicit expression for Dyson equation can then 
be obtained by taking the symmetric sum of Eqs. ([20]) 
and (jA6|) . and afterward, by expanding the exponential 
function and using the binomial theorem: 

xa^9^"-^«,„(R)ai"-^9^5f„'^(R). (A7) 

Note that the difference of Eqs. and yields 

another equation which may be useful in solving Eq. (|A7p 
(e.g., in the case of a quadratic potential, see Sec. 

+ 00 /,4 \ " 2ri+l , N^„+„ 



, 4 / ^ p!(2n+ 1 -n)! 



xdldl-+'-^Vrr.{^)df+'-^dl.~g^^{^). (A8) 



APPENDIX B: MODIFIED VORTEX WAVE 
FUNCTIONS 

Our aim in this appendix is to prove expression 
Let us analyze first the following differential operator 



+ 00 . ,2 ^ p 



p=0 ' 



(Bl) 



Applying this to a function /(R) and introducing the 
Fourier transform of /, we get 



±^ 1 / /? 



p=0 ' ^ 



(B2) 



/0E,^(^O /»e--(B3) 
J-^me^^-^e^'^^'/\ (B4) 



(2^) 

Using the inverse Fourier transform, we have 

"^^^ „iq(R-u 



0[/(R)]= /d^u/(u) 



27r)2 



g„q-(R-u) plW/i_ 



(B5) 



The integral over q is formally divergent. We circumvent 
this problem by introducing for a while the parameter 
^ = — /^/4 > 0. The calculation of the resulting Gaussian 
integral can then be easily done, which finally yields 

6[/(R)] = -|^/(u)e(^-")^/'l. (B6) 

We deduce from this calculation that the operator O is 
nothing but a convolution operator with a Gaussian ker- 
nel. We now apply it to /(R) = ^*^ j^{r')^,„,.Rir). Using 
formula (jBOp and the explicit expression ^ of the vortex 
wave functions, we have 



0[*:„,R(r')*™,R(r)] = - 



,-(R-c)ViBgi[(2R-c)xd]-J/i^ 

27r/|TO! 



fd^V 




m 









Jr,~2[R-c]f /2ll 



(B7) 



where we have done the change in variable 77 = u — c + 
id X z with d = (r' — r)/2 and c = (r' + r)/2. We are 
again in presence of a formally divergent integral. As just 
above we introduce the parameter ^ and use the following 
trick to perform the Gaussian integral over -q in Eq. (|B7[) : 



2/2 



,(r,-2[R-c])V2i| _ 









m 

-(r,-2[R-c])V8C 






[ 8C J 






g,n 


fd^V 


-S'VVSC „-(»7-2[R-c])V8? 












s- 



(B8) 



The remaining Gaussian integral (jB8[) can now be 
straightforwardly evaluated (note that the contours of 
integration can be deformed to the real axes using the 
analyticity property of the integrand). We finally find 



0[*:.,R(r')*™,R(r)] 



-(R-c)^/i|jgi[(2R-c)xd]-z/;^ 



ds 



1 + s 



(B9) 



s=0 



Inserting the definitions of the parameters c and d in 
terms of the positions r and r' into Eq. (|B9p . we directly 
arrive at expressions and 



APPENDIX C: CHECKING THE VORTEX 
FORMALISM: CASE OF A ID PARABOLIC 
CONFINING POTENTIAL 

1. Standard derivation 

In the particular case of a ID parabolic potential given 



by 



2^2 



V{x) — -m*u;QX 



(CI) 
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the wave functions and the energy spectrum of Hamilto- 
nian ([T]) can be found by solving directly the Schrodinger 
equation using well-known standard methods. The rele- 
vant quantum numbers appear to be a positive integer n 
which labels the Landau levels, and a continuous quan- 
tum number py playing the role of momentum in the y 
direction. In the Landau gauge A = Bxy, wave functions 
and energy spectrum read, respectively. 



^-[x-{LjJn)pyL^f/2L^ 

*«P,(r) = , = — Hn 

\/2"+in! 7r3/2L 

ipy 



L 



npy 



xe 



V{pyL^ 



(C2) 
(C3) 



where f2 = 



' Wq and L 



are the renor- 



malized cyclotron pulsation and magnetic length, respec- 
tively, and H„ denotes the nth Hermite polynomial. 

In absence of Landau level mixing, one has to con- 
sider Wc ^ Wo, keeping all terms of order ujq/lUc, and 
neglecting higher powers of this ratio. Thus we have 
n « uJc+uJo/2uJc and L^Ib- From Eqs. (IC2| and ([CS]) . 
the Green's function thus reads in this limit of negligible 
Landau level mixing 



-t-oo 



G«'^(r,r',a;) = ^ dpy 



xHn 



2"+in!7r3/2/B 

it' — n. 

Hr, 



X Pyl%\ „ fx' Pylg 



h 



(C4) 



with 



Enp ~ tuJc ( + 2) ^ ^^°2tj" ^ 2 ^ ^ ^(Pv^b)- 



(C5) 



where c and d are defined in Appendix |B1 Eq. (jC6p is 
rewritten as 

+00 



dX 



1 



, 2^1^ uj — E„ 

rn=0 •' o " 



Vm{X) ± i5 



ds 



=e B e 



' 



(C8) 



s=0 



It can be checked that the following algebraic relation 
holds: 



9s™ 



1 



zc Be 



e B H, 



J ,s=0 

x' -X 



Inserting formula (jC9p into Eq. (|C8[) and reintroducing 
the variables r and r' everywhere in place of c and d, we 
find that the Green's function finally reads 



+00 „ 



1 



dX 

1% - E„i - Vjn{X) ± iS 



xe 

x-X 



B 



Hr, 



2™+im!7r3/2ZB 

x'-A^ ,sv-^ 



(CIO) 



Introducing py — X/lg and expliciting the term Vm{X) 
by inserting expression (jCip into definition (fT8|l of Vm, 
we see that expression ()C10P corresponds exactly to Eq. 
(|C4[) up to a phase factor exp [i{xy — x'y')/2Pg\ which 
comes from the fact that we work here within the vortex 
formalism in the symmetric gauge, and not in the Landau 
gauge. 



2. Derivation within the vortex formahsm 

Now, we show how one can recover the Green's func- 
tion (|C4p of a ID parabolic potential from the vortex 
formalism. We start with expressions pT|) and ([24|) and 
exploit the fact that the effective potential ■;;,„ is inde- 
pendent of the variable Y 



G 



R,A 



+00 „ 

(r,r',c.)=^ /. 

n "J 



dX 



1 



Tn=0 B 



X dYe- 



Em - Wm(A) ± i6 

[^;.,R(r')*™,R(r)] . (C6) 



The integral over Y can then be performed exactly mak- 
ing use of expressions (|22|) and ([23|) . Considering that 



IdYe 




"(C7) 



APPENDIX D: SOLVING THE DYNAMICAL 
EQUATION FOR POTENTIAL LINES 

Differential Eq. pTjl is second order in the derivative 
with respect to E, but first order in E. It will obviously 
become second order in r and first order in the deriva- 
tive with respect to t by going to the Fourier component 
F^'^{t). So, in order to solve Eq. (pij) . we write 



-iEt 



(Dl) 



f^,'^(E)^ JdrF„Y{r)e 
and substitute this form into Eq. (PT|) to get 



xe-'^^ = 



(r) 
1. 



ujm ± i(5 — i^T — TjT^ — il'^T'^ + l)";— 

dr 



(D2) 
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Doing an integration by parts, we have 



1= [-»(l + 7r^)F^'^(r), 



iET-\+°° 



-iEr 



UJm as + i^T — 77T^ + i (l + 7T^) 



dT 



dr e 



Finally, taking the Fourier transform of this equation, we 
find that F^'^{t) is governed by the first-order differen- 
tial equation 



i [l -I- 7t2] 



dF^-^^{r) 



dr 



+ [ura - vr^ + nr ± iS] F^^'^{t) ^ S{t), (D4) 
provided that the integrated term in Eq. (jD3p vanishes. 



of the boundary term at infinity. For 7 > 0, it can be 
easily checked that the function 



f^-^{E) = T^ Jdt 



cos(^t) 



xe 



-i(E+7]/t)r(t) 



(D9) 



(DIO) 



with the function r(t) given by 

r(i) = tan {y^t) , 

is a solution of Eq. (|3ip . Here integral (|D9[) is defined 
in the sense of Cauchy principal value for the points 
y/jt = 7r/2 -I- riTT. This provides the exact result ([35]) 
for the vortex Green's function of an arbitrary quadratic 
potential in the lUc 00 limit. 



{l + jT')F,f^{T)^0 (D5) 
when T ztoo. Equation (|D4p is readily solved by 



V^l + it' 



: exp 



T] 

i {u},n + ?7/7 ± "iS) t{T) - i-T 

7 

(D6) 



where we have introduced 
1 



tir) 



V7 



arctan {^/JT) . 



(D7) 



Here 6{t) is the Heaviside function. For 7 < 0, one must 
understand that 



1 



arctan {y/jr) 




defined for \/— 7|t| < 1. The variable t in expression (|D6P 
plays actually the role of the time since it is conjugated 
to the energy lu which enters into the expression via the 



quantity Cbm, see definition 



Because the solutions 



of the homogeneous equation do not respect the time 
causality, one has only considered the particular solution 
of inhomogeneous equation (|D4[) . 

For 7 < 0, solution (|D6p fulfils requirement (|D5P for 
any value of the parameter 77 (for the case 7 = 0, con- 
dition (|D5[) is obeyed with the help of the infinitesimal 
quantity ztiS, while for 7 < we have F^-^ = for 
7|t| > 1). However, for 7 > 0, we note that con- 
dition (jPSp is not satisfied, so that expression (|Dip to- 
gether with formula (|D6p does not yield a solution of the 
initial Eq. (PT|) . Nevertheless, the solution of Eq. ([31]) 
for 7 > can be inferred from result (|D6[) by noting that 
the problem actually originates from the saturation of 
the function t(T) when r — > ±00. Indeed, by consider- 
ing t instead of r as being the relevant variable and by 
extending its domain of definition to the whole real axis, 
we can exploit the infinitesimal quantity ±iS to get rid 



APPENDIX E: CHECKING THE VORTEX 
FORMALISM: CASE OF A 2D PARABOLIC 
CONFINING POTENTIAL 

Recovering the set of two discrete quantum numbers 
for the circular confinement potential from the use of a 
basis of states which is characterized by both discrete and 
continuous quantum numbers appears in principle to be 
a very challenging task. We show that the quantization 
of the confining potential appears in the vortex Green's 
function formalism through a rather different way from 
the usual derivation in the wave function formalism. 



1. Standard derivation 

To benchmark our results for the Green's hmctions, we 
shall compare the general expression derived in Sec. IIVI 
from the use of the vortex states formalism with the exact 
solution for a circular confining potential. The potential 
profile given by 



(El) 



leads in a homogeneous magnetic field to the well-known 
Fock-Darwin spectrum 



Eni = n 



N + 1 



huJr 

2 



(E2) 



where n 
0,±1,±2, 



0,1,2,... is a positive integer and I 
a positive or negative integer. Here 



+ AijJq is the renormalized cyclotron pulsation. The 
normalized wave functions associated with the energy 
spectrum (|E2[) are written in polar coordinates r = (r, 6) , 



1'n,i(r) 
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where Ln{z) corresponds to the gcnerahzed Laguerre 
polynomial of degree n, and L ~ \J h/nVVl is the renor- 
malized magnetic length. 

The local density can be directly calculated from the 
knowledge of the energy spectrum and the exact wave 
functions, and is given by 

+OC +00 

"W = E E "F(£;„OI*«Kr)r. (E4) 

n— — 00 

The method of projection onto a given Landau level is 
again obtained by considering S> wq, keeping terms of 
order loqIloc- This is equivalent to taking ^ 00 with 
Ib finite. We thus have O « Wc + Swq/wc and L Ib, so 
that the energy spectrum becomes 

/ 1 \ 2^0 / l + , , 

Eni sa hujc ( ™ + 2 ) ^ lJ~ \ 2~ ) ' 

with m = [n + \ — I) / 2] > the Landau level index. Ac- 
cording to the second term in the right-hand side of Eq. 
(IE5p ; the Landau levels are generally nondcgenerate as a 



result of the circular confining potential characterized by 
the frequency coq tUc- 

If we restrict ourselves to the lowest Landau level con- 
tribution to the local density for the sake of simplicity 
and consider the absence of Landau level mixing, the ex- 
act local density gets simplified into 

xe"^. (E6) 



2. Derivation within the vortex formalism 

The different parameters for the circular confining po- 
tential are ( = /^m*cjQ = HllIq/uJc, 7 = and 
T]{r) = lg( I VV(r)|^ /8. Using these values and the gen- 
eral formula for the local density pS)) obtained from the 
vortex formalism, we get for the lowest Landau level con- 
tribution (m = 0) to the local density 



n(r) = 



1 



2ttI% 
1 



2ttI% 



Im / dt 



T 



sinh [nTt 
T 



- + Im dt , , 
2 Jo sinh [nTt] 



exp 



cxp 



2 

hwc 



+ 00 



1=0 ^ ^ 



00 rp 

dt—- 

sinh [TrTtl 



^0 



^+2P ' 



E n \2ll 

1=0 ^ ^ 



-ilthuJg/LJc 2l| 



sm 



,---(/ + l)-^)t 



(E7) 
(E8) 

(E9) 



Using the integral 



-f-oo 

dx 



sin(aa;) tt 



/o sinh(6x) 
local density (|E9p is rewritten as 



Finally, by noting that the first term in the right hand 
side of Eq. (|E11|) can be written as 



1 1 

2^n \2il 

1=0 ^ ^ 



(E12) 



7i(r) 



1 



1 IP^l f 



2nl% [2 + 2E;! [211 



X tanh 



,- — -il + l)-J>]t 



(Ell) 



we arrive at formula (|E6|) for the local density. This 
establishes the exact equivalence of general formula ((48|) 
and of Eq. (jE6| in the particular case of a circularly 
symmetric confining potential. 
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